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Abstract 

We represent the slow, glassy equilibrium dynamics of a line in a two- 
dimensional random potential landscape as driven by an array of asymp- 
totically independent two-state systems, or loops, fluctuating on all length 
scales. The assumption of independence enables a fairly complete analytic 
description. We obtain good agreement with Monte Carlo simulations when 
the free energy barriers separating the two sides of a loop of size L are drawn 
from a distribution whose width and mean scale as L}^^ , in agreement with 

recent results for scaling of such barriers. 
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Slow dynamics is perhaps the most significant characteristic of the glassy state of matter, 
affecting essentially all experimental measurements. An intuitively appealing picture, which 
explains this remarkable slowing down, is that the configuration space of a glass consists of 
many nearly degenerate free energy minima, separated by high potential barriers The 
dynamics is dominated by transitions between configurations whose free energy difference 
/\E is less than the thermal energy A;_bT, and which are separated by free energy barriers 
Eb ^ ksT. In a genuine glass such degeneracies occur on all length scales L. The de- 
pendence of AE and Eb-, or more precisely, of their probability distributions, on the linear 
extent L has become the focus of great theoretical interest . If transitions mainly occur 
between pairs of low-energy configurations which can be regarded asymptotically indepen- 
dent for large L, the simple model of a gas of fluctuating two-level systems allows for a 
fairly complete description of the dynamics. In general, however, a much more complicated 
hierarchical interdependence of transitions on different length scales may take place 

The subject of this Letter is an elastic line (henceforth called an interface) in a two- 
dimensional random potential landscape This system combines a fair degree of 
realism e.g. as a model for a domain boundary or a magnetic flux line trapped between two 
copper-oxide planes in a dirty high-temperature superconductor |^,^ , with a simplicity that 
has allowed a substantial body of knowledge to accumulate over the past decade [|TOi|TT| . An 
almost degenerate two-level system in this case is simply a loop: a segment of the line between 
two points, which can flip between two low energy paths (valleys in the potential landscape) 
separated by a barrier (a mountain in the landscape), this is illustrated in the inset of Fig. |l|. 
It has been well established 0,§J3] ^^^^ transverse size of such a loop scales as A/i oc L'^ 
with = 2/3, while the free energy differences between the two valley configurations, AE 
are distributed with mean zero and variance {AE"^) oc L^^, where 6 = 1/3. Recent work ^ 
provides evidence that the barriers between such configurations are distributed with mean 
{Eb{L)) oc and the variance {E"^) oc L^^ with likely logarithmic corrections. 

These results, when combined, provide all necessary ingredients for developing the dy- 
namic description of the line under the assumption of asymptotic independence of the transi- 
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tions within large loops. In this Letter we outline such a description, and show that it agrees 
well with dynamical Monte Carlo simulations. We first numerically confirm that nearly de- 
generate paths form loops of various sizes, and at first stage we analyze the dynamics of 
loops of a fixed length. The dynamics of a single loop can be described as a sequence of 
flips where the interface position moves from one arm of the loop to another. By studying 
loops of different lengths we determine how the fiipping-rate distribution depends on the 
size of the loop. At the second stage we study the fluctuations of a low-energy interface. 
Using numerical simulations we determine the time-dependent fluctuations of the interface 
position, which we compare with an analytic model that assumes that interface fluctuations 
are due to flipping loops. Since parameters describing loops are determined at the first stage, 
there are no adjustable parameters left at this point. We find that the prediction of the loop 
model is in good agreement with the simulations, supporting the conjecture that interface 
dynamics is due to fluctuating loops. 

The numerical model: We study a model that excludes overhangs so that the interface 
height h{x) is at all times a single-valued function of the spatial coordinate x. For the 
numerics we use a lattice model where the interface is discretized to have a unit slope between 
the lattice points, \h{x) — h{x+l) \ = 1, and use fixed boundary conditions /i(0) = /i(-Lo) = 0, 
where Lq is the interface length. The Hamiltonian of the lattice model for a particular 
realization of the random medium is 

H,[h] = Y^fxih{x),x), (1) 

x=l 

where h{x),x = 1, . . . ,Lo represents the interface, and /i is the potential landscape. The 
random potential is uncorrelated and uniformly distributed over the range < fi{h,x) < 1 
1^ . This Hamiltonian belongs to the universality class of a directed polymer in a random 



medium (DPRM) [0,|lH], and in this context the requirement of unit slope between lattice 



points effectively adds an implicit line tension [|TT| . The dynamics is implemented in a spirit 



similar to that of previously introduced mappings onto spin chains |T^,|TB[; the dynamics is 
described by the master equation based on the transition probability P[h{x) h'{x)] = 



exp{j3{Hfj,[h] — Hfj,[h'])/2}dt, where j3 is the inverse temperature, dt an infinitesimal time 
interval, and h{x) and h'{x) are two interface configurations that differ in only one position. 
Numerically the dynamics of the master equation is exactly modeled by an algorithm that 
uses time steps sampled from a Poisson distribution [|T6|. We have used /3 = 2 in the 
simulations that are presented in this Letter. In general the computation time grows very 
rapidly with increased f3. 

Single loop statistics: In order to explore the ideas of loop-based dynamics we use 
the free energy landscape to define loops. The free energy of a point (x, h) is defined as 
F{x, h) = —(3~^ ln(P(a;, h)), where h) is the probability that an interface with fixed ends 
{{x,h) = (0,0) and (Lo,0)) crosses this point. The probability P{x,h) is easily calculated 



using transfer matrix methods ||TT|. A lattice point is defined to be part of an island if it 
is not part of any interface that includes only points with P{x,h) > 0.1. The level 0.1 is 
chosen to obtain well defined islands, but the specific value does not greatly affect the final 
results of the analysis. Interface segments encircling an island form a loop, and we measure 
the loop size in terms of its length L = Xr — xi + 1, where {xr, h^) and (x;, hi) are the right 
and left ends of the island, respectively. The center height of the island, his, is defined as 
the /i-component of its center of mass, where every lattice point in the island is regarded as 
a mass point. We say that the loop surrounding the island changes its state (fiips) when the 
interface height, averaged over the island length, crosses his. This definition is convenient 
although it sometimes catches events that would not intuitively be considered as fiips. 

In order to study the relationship between static and dynamic scaling, and to unam- 
biguously determine the numerical values of the amplitudes to be later used in the fit, we 
collect statistics of the dynamics of individual loops. A loop is inscribed in a bounding box 
which sets the boundaries within which the interface is free to move (showed with dashed 
lines in Fig. |1|). The interface is constrained to pass through the left and right corners of 
the bounding box, (x/ — 2, hi), and (x^ -l- 2,hr). We collect statistics of the time between 
consecutive fiips of the loop. In Fig. |T] we plot a typical example of the measured probability 
that a loop stays in the upper (lower) state for at least time t. The decay of this probability 



is exponential for sufficiently large t, which is consistent with the two-state model discussed 
below (Eq. (Q)). The deviation for small t is probably due to our simplified criterion for 
flipping the loop. Before making a least-square fit to an exponential decay we exclude the 
part of the short-time data that strongly deviates from the expected form. We also exclude 
the largest times (0.1% of the data) since the statistics of these very rare events is poor. 

From this we obtain the characteristic decay times and r_+ for flipping the loop from 

its upper to lower state and vice versa. 

For an individual loop we collect data from 10,000 flips, and make the fit described 
above. The rate constant, F, characteristic for the loop, is calculated by F = r^i -|- rr+ 
(we assume that the free energy difference between the arms of the loop is unimportant). 
By collecting statistics from 1000 loops of the same size L we find that F is log-normally 
distributed. This is consistent with a simple activated behavior, F = Fe"'^'^, where A is a 
Gaussian-distributed energy barrier separating the two sides of the loop, and F is a constant 
setting the unit of time. In Fig. ^ we show fits to the log-normal distribution by plotting 
Q^^{P(r)) against ln(F), where Q is the complement of the cumulative normal distribution, 
and -P(F) is the measured cumulative probability of the rate constant F (hence, a straight 
line would correspond to a log-normal distribution). The distribution is characterized by 
the average barrier height A(L) and the standard deviation (j{L). The values for (3a{L) are 
obtained as the standard deviations of ln(F), which also give the inverse slopes of the lines 
in Fig. H. Similarly, the crossings between the lines and the zero axis represent the averages 
of ln(F), which give values for ln(F) — /3A(L). The average barrier height A(L) and the 
standard deviation <j{L) scale with loop size as L^/^ ||^. This allows us to determine the 
time scale F and the multiplicative constants in A(L) and cr(L). The mean and standard 
deviation of ln(F) are plotted as functions of L^/^ in Fig. ^ conffiming the L^/^ scaling, and 
giving for /5 = 2 the parameter values ln(r) = 2.2, /3A(L) = 3.0L^/^ and (3a{L) = 0.28L^/l 

Full interface numerics: The same numerical algorithm is used for studying the dy- 
namics of the full interface but keeping only the ends of the interface fixed. We start the 
interface from a random initial state (chosen from the equilibrium ensemble), measure how 
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the height hc{t) of the center point of the interface varies with time, and collect statistics of 
6h'^{t) = [hc{t) — /ic(0)]^. The equilibrium ensemble is used to normalize 6h'^{t) with respect 
to 5/i^(oo). The data points in Fig. ^ show the results of simulating interfaces of lengths 
20,40,60, and 128 for /3 = 2, where 6h'^{t) has been averaged over 20,000 realizations of the 
random medium. 

The analytic model: We model the loop dynamics using a two-state model, where the 
probability of a given loop being in the upper or lower state, respectively, is given by P+{t) 
and P-(t). The time development of these probabilities is governed by the coupled differen- 
tial equations 

(2) 



dt 



r+_p+(t)-r_+p_(t). 



The fluctuation in the height of the center position of a loop of width wLF'^^ is given in the 
two-state model by 

5hl,{t) = {[Kit) - K{0)f) = ^L^/3(l - e-^*), (3) 

where we assumed for simplicity r_|__ = r__|_ = T/2, i.e. that the two arms of the loop are 
exactly degenerate. 

We consider an interface of length Lq, and denote the average barrier height of the largest 
loops (i.e. loops of size Lq) by Aq, and its standard deviation by ctq- We assume that the 
number of loops of a given size scales as L"^, and that the fluctuations of different loops are 
independent and additive. We find that the total fluctuation of the interface, as implied by 
loop dynamics, is given by 



SK{t) dy 



1 ,,2 



=e 2' 



5/1^(00) J-00 v27r "'0 

The first integral corresponds to integrating over loops of fixed length but variable rate 
constants, and the second integral adds up the contributions of loops of different lengths. 
The infinite time fluctuations 6h'^{oo) are given by the equilibrium result, and scale as Lq^^ 
i- 



The two-state model (Eq. (P) implies that the natural time scale in the problem is 
given by t = (r)-i ln(2) exp(/?Ao/2i/4) in the sense that at time t = t the fluctuations have 
reached approximately 50% of the equilibrium value; more precisely 0.5 < ^1^2^-^ ^ 0.58 for 
all /3, Aq, and ctq. Hence, at low temperatures interface dynamics are exponentially slow as 
is typical for glassy systems [jl^ . 

The solid lines in Fig. ^ are the results of the two-state model using the parameters 
determined by studying individual loops. The time scale F determines the position of the 
curves, the average barrier height, Aq, affects both the position and the slope of the curves, 
and the standard deviation of the barrier heights, o"o, has only a minor effect on the results 
for small (To/Aq (taking the limit (To/Aq would not significantly change the results). 
Considering that there are no adjustable parameters, the agreement with the results of 
the numerical simulations in Fig. ^ is quite good. This suggests that the hypothesis that 
interface dynamics is due to fluctuating loops is indeed valid. The deviations in the short 
time behavior are, we believe, due in part to more complicated loops on small length scales, 
which cannot be described as independent simple loops. Another contribution to the short 
time deviations is that in the lattice model the arms of a loop can have a non-zero width (i.e. 
two adjacent interface positions are not separated by a barrier), which naturally enhances 
fluctuations in short time scales. The interpretation that small time differences are due to 
discretization and finite size corrections is further supported by the fact that the deviations 
are smaller for larger system sizes. Another possible source of deviations is logarithmic 
corrections to the scaling forms of A(L) and cr(L), however, our data on individual loops is 
insufficient to determine these corrections. 

In conclusion, we have studied the dynamics of a one-dimensional interface in a two- 
dimensional random medium. We have shown that the dynamics can be understood quanti- 
tatively in terms of loops formed by nearly degenerate interface paths. At low temperatures 
the dynamics is exponentially slow, which is typical for glassy systems p|,|^,p^ . 



7 



REFERENCES 

(a) P];gggnt address: Prudential Home Mortgage Company, 8000 Maryland Ave., Ste. 1400, 
Clayton, MO 63105. 

[1] M. Mezard, G. Parisi, and M. A. Virasoro, Spin Glass Theory and Beyond, World 
Scientific, Singapore (1987). 



[2] 
[3] 
[4] 
[5] 
[6] 
[7] 



[9 
[10 



[11 
[12 

[13 
[14 
[15, 



D. S. Fisher and D. A. Huse, Phys. Rev. B 38, 386 (1988); 38, 373 (1988). 

D. A. Huse and C. L. Henley, Phys. Rev. Lett. 54, 2708 (1985). 

D. A. Huse, C. L. Henley, and D. S. Fisher, Phys. Rev. Lett. 55, 2924 (1985). 

L. loffe and V. M. Vinokur, J. Phys. C 20, 6149 (1987). 

D. S. Fisher, M. P. A. Fisher, and D. A. Huse, Phys. Rev. B 43, 130 (1991). 

D. S. Fisher and D. A. Huse, Phys. Rev. B 43, 10728 (1991). 

L. V. Mikheev, B. Drossel, and M. Kardar, Phys. Rev. Lett. 75, 1170 (1995). 

P. W. Anderson, B. I. Halperin, and C. M. Varma, Philos. Mag. 25, 1 (1972). 

M. Kardar, in Les Houches 1994, Session LXII, fluctuating geometries in statisti- 
cal mechanics and field theory, edited by F. David, P. Ginsparg, and J. Zinn- Justin 
( Ihttp: / /xxx.lanl.gov/lh94/e-book/| or |cond-mat/9411022|) . 



T. Halpin-Healy and Y. C. Zhang, Phys. Rep. 254, 215 (1995). 

A Gaussian distribution with the same mean and standard deviation gives results that 
marginally differ from those of the uniform distribution. 

H. Yoshino, unpublished ( |cond-mat / 95100"2^ ) . 

L.-H. Gwa and H. Spohn, Phys. Rev. A 46, 844 (1992). 

H. C. Fogedby, A. B. Eriksson, and L. V. Mikheev, Phys. Rev. Lett. 75, 1883 (1995). 



[16] K. Binder, Sec. 1 in Monte Carlo Methods in Statistical Physics, edited by K. Binder 
( Springer- Verlag, Berlin 1979). 

[17] K. H. Fisher and J. A. Hertz, Spin Glasses (Cambridge University Press, Cambridge, 
1991). 



9 



FIGURES 




rescaled time 

FIG. 1. Typical time dependence of the probability that a loop stays in the upper (lower) state 
for at least a time t. The data is based on 10,000 crossings in each direction for an island of size 
L = 10. Note that the time scales are normalized to make the right-end value 0.1% of the initial 
value. The inset schematically shows a loop formed by two nearly degenerate alternative paths of 
the interface. 




FIG. 2. The cumulative distribution, P(ln(r)), and the cumulative normal distribution Q are 
used to illustrate the fit to a log-normal distribution of the T data (solid lines) for the loop sizes 
L = 1, 2, 4, 8, 12, 16, 20 (from right to left). The statistical data is based on 1000 loops of each size. 
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FIG. 3. The parameters A(L), cr{L), and T are determined by fitting (3A{L) — ln(r) to the 
mean, and (3(t{L) to the standard deviation of the ln(r) data. The error bars show the 95% 
confidence interval of the statistical error. 
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FIG. 4. The points show the time dependence of 5h?{t) for interfaces of lengths 20, 40, 60, and 
128 averaged over 20,000 samples each. The solid lines show the result of the loop model (Eq. ^) 
where the parameters of the single-island analysis have been used. 
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